Mon. Not. R. Astron. Soc. 000, [JHS] (2009) Printed 7 May 2010 (MN WI^ style file v2.2) 



Cosmic Rays from Pulsars and Magnetars 



o 



O 



> 

m 

o 
o 



o 
o 



X 



Jeremy S. Heyl^^, Ramandeep Gill^ and Lars Hernquist^ 

^Department of Physics and Astronomy, University of British Columbia, Vancouver, British Columbia, Canada, V6T IZl 
'^Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge MA 02138 
^ Email: heyl&phas.ubc.ca; Canada Research Chair 



Accepted — . Received — ; in original form 



ABSTRACT 

We compare the expected abundance of cosmic-ray electrons and positrons from pul- 
sars and magnetars. We assume that the distribution of infant pulsars and magnetars 
follows that of high-mass stars in the Milky Way and that the production rate of cos- 
mic rays is proportional to the spin-down and magnetic-decay power of pulsars and 
magnetars, respectively. In combination with primary and secondary cosmic-ray lep- 
tons from other sources (especially supernova remnants) , we find that both magnetars 
and pulsars can easily account for the observed cosmic-ray spectrum, in particular 
the dip seen by HESS at several TeV and the increase in positron fraction found by 
PAMELA. 
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1 INTRODUCTION 

Recent measurements of cosmic rays from en ergies of one 
GeV to one TeV have sho wn a spectrum (jChang et al.l 
120081: Aharonian et al.l l200g ) and abundance of positrons 
( Adriani et al.l [2009|) that are difficult to reconcile with 
astrophysical sources. These observations have naturally 
stirred up great interest among the particle physics com- 
munity and e specially people who study dark matter (e.g. 
lHoopeiil2009l ). Essentially the problem is that the observed 
spectrum is much harder than that expected to be produced 
by astrophysical sources such as supernova shocks or pulsars, 
so it is quite natural to look at more exotic possibilities, in- 
cluding decaying dark matter. On the other hand, several 
works have also examined the question of whether or not 
astrophysical sources can account for the observed spectrum 
and positron abundance by including a reasonable model for 
the diffusion of cosmic rays through the Galaxy. In this Let- 
ter we consider the expected abundance of cosmic-ray elec- 
trons and positrons (hereafter referred to as cosmic rays) 
produced by magnetars and pulsars. 

Magnetars (IDuncan fc Thompson! 1 19921 : 

iThompson fc DuncanI 1 19951 ) are ultramagnetized neu- 
tron stars with fields B. ~ lO " - 10^*^ G, fueled 
by e i ther magnetic field decay (jThompson fc DuncanI 



1 19961: Hevl fc Kulkarnil 19981 ) or residual thermal energy 



( Hevl fc Hernguistl Il997al lbl). Magnetars are thought 



to 
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produce co pious numbers of electron-positron pairs both 

in bu rsts 

2005bl) 
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Beloborodov 



Hevl fc Hernguistl 

and in quiescence (JHevl fc Hernguistl l2005al : 

ThompsonI |2007| V Moreover, the local 



radiation from these pairs is thought to power the bursts 



from soft-gamma repeaters (JThompson fc DuncanI Il995l ) 
as well as the recently discovered non-thermal emission 
from soft-gam ma repeaters and anomalous x-ray pulsar s 
in quiescence (|Mereghetti et all [20051 : iKuiper et al.ll2006l ). 
Consequently, magnetars seem a natural candidate to make 
an important contribution to the observed abundance of 
elec trons and positr ons observed at Earth. 

IProfumd (|2008l ) presented a detailed calculation of the 
expected spectrum of cosmic rays in various models for 
their production from pulsars. His calculation focused on the 
population of gamma-ray pulsars. The situation for magne- 
tars is somewhat different because these objects are much 
rarer than gamma-ray pulsars; the birth rate for magne- 
tars is roughly ten times smaller than that of pulsars. Con- 
sequently, the cosmic ray flux that we observe can come 
both from the magnetars that are active today as well as 
from more ordinary looking objects that may have been 
active magnetars in the past such as RXJ 0720.4-3125 
iHevl fc Hernguistlll998l : iHevl fc Kulkarnil 1 19981 ). To calcu- 
late the expected abundances of cosmic rays from magne- 
tars, we generate a Monte Carlo simulation of the magnetar 
population over the past ten million years and compare it 
with a similar simulation of the pulsar population over the 
past million years. 



2 THE MODELS 

2.1 Emission 

To model the cosmic ray emission from pulsars we assume a 
production rate always proportional to the spin-down powe r 
by magnetic dipole emission. In contrast to lProfumd (|2008r ). 
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we do not make any distinction between objects that we ex- 
pect to be gamma-ray pulsars versus the rest of the pulsar 
population. To specify the model, we assume that all the pul- 
sars are born with a period of 49 milliseconds and spin down 
to one second by the time we observe them. For the magne- 
tars, we take the cosmic ray emission to be proportional to 
the magnetic field deca y power. We use the calculations of 
iHevl fc Kulkarnil (|l998r ) to follow the evolution of the mag- 
netic field starting at 10^® G at the pole on the surface. We 
assume that only the energy in the external magnetic field 
is converted into cosmic rays and that the external field is 
strictly dipolar (if the internal field is uniform, there is a 
factor of fifty more energy within the star). For ease of com- 
parison, the initial spin period of the pulsar population is 
chosen such that the total spin-down energy in the pulsars 
equals the total magnetic field energy of the magnetarslj The 
magnetars and pulsars may inject high-energy electrons and 
positrons into interstellar space. We assume that the fraction 
of the magnetic field or spin-down energy converted to pairs 
is given by e„i and ep respectively. These are determined 
by insisting that total electron-positron cos mic-ray abun- 
dance from other expected sources (Dclahav c et al.l [2010|) 
when combined with either the magnetars or pulsars does 
not exceed that observed by Fermi around 300 GeV. 

An important difference between magnetars and pulsars 
is the evolution of cosmic ray generation with time. In the 
case of a pulsar the bulk of the cosmic rays are made during 
the very early life of the pulsar as it is spinning down most 
rapidly. On the other hand, we assume that magnetars gen- 
erate cosmic rays from magnetic field decay, which is also 
thought to power thermal and nonthermal photon emission. 
The magnetic field decay timescale is typically several thou- 
sand years so the emission of cosmic rays from the magnetar 
is of longer duration. This difference shows up in the ob- 
served spectrum of cosmic rays from a nearby magnetar as 
compared to a nearby young pulsar (Fig. [T]). In the case of a 
pulsar, one sees a relatively smooth spectrum up to the en- 
ergy of the electrons or positrons that have had time to cool 
since the birth of the pulsar; at this e nergy there is a rela- 
tively sharp cut off. In the case of the iProfumd ((2003) cal- 
culation this cut off is infinitely sharp. A magnetar releases 
its magnetic energy gradually, so the cosmic rays that we 
observe at Earth have a variety of ages, and the energy cut 
off is much less sharp. However, this is only a second-order 
effect in determining the shape of the spectrum from an 
ensemble of objects; the largest differences in the expected 
spectrum from a collection of magnetars versus pulsars owes 
to the relative ages of these two populations. 

A second important issue is the spectrum of the cos- 
mic rays where they are produced. In both cases we as- 
sume that the number of cosmic ray particles decreases as 
7~'^'^ (7 is the Lorentz factor of the cosmic ray), as sug- 
gested by observations of synchrotron emission from shocks 



^ We note that these choices, while somewhat arbitrary, are mo- 
tivated by observed characteristic properties of young and old 
pulsars and magnetars. Moreover, variations in pulsar spin pe- 
riods and the magnetic energy available to produce cosmic rays 
from magnetars affect only the normalization of the cosmic ray 
flux seen at Earth, not their observed energy spectrum, which is 
mostly sensitive to the distribution of sources in space and with 
time and the spectral index of the intrinsic emission process. 



(jSalmonson et al.ll2006f ). Furthermore, several strongly mag- 
netized neutron stars appear to have emission similar to 
pulsar wind ne bulae (^Rca ct al.l l2009l : iHelfand et al.l l2003l : 
iGonzalez fc Safi-Harb 2003 ). The apparent absence of such 
structures around the classical magnetars (the soft-gamma 
repeaters and anomalous x-ray pulsars) may simply owe to 
the difficulty of detecting them at their relatively larger dis- 
tances. Perhaps these nebulae are powered by pulsar spin- 
down; however, the efficiency of the conversion of spin-down 
energy to nebular emission appears to be higher than in 
typical puls ar wind nebulae , so magnetic field decay may 
play a role (JRea et al.ll2009l ). Moreover, magnetars may be 
born w ith spin rates characteristic o f typical pulsars (or even 
faster; iThompson fc DuncanI 119931 ) so in their youth they 
may emit pairs in a similar manner as pulsars (of course, the 
time evolution of the emission would differ but this only has 
a second-order effect on the observed spectra). Therefore, 
with supporting hints and without definite observational ev- 
idence to the contrary, we assume that the magnetars emit 
a similar spectr um to pul s ars. T his spectrum is softer than 
that adopted bv lProfumd (|2008l ) who used oc 7"^'^ through 
7~'^. The spectrum of particles produced is taken to extend 
from 7 = 1 to 7 = 10^. 

2.2 Diffusion 



Like [ Profumd (|200a ). we use the model of lAtovan et al.l 
(|l99a) to calculate the diffusion of cosmic rays through 
the Galaxy. We will briefiy summarize previous re sults and 
the di fferences between our treatment and that of lProfumd 
(|2008h . 

Because of the Galactic magnetic field, cosmic rays can- 
not travel directly from their source to Earth. Rather, they 
diffuse through the Galaxy with higher energy particles hav- 
ing larger Larmor radii diffusing more quickly than lower 
energy particles. Meanwhile, as the particles gyrate about 
magnetic field lines, they also lose energy. If the region over 
which the particles diffuse is larger than the volume con- 
taining the important sources, it is appropriate to assume 
that the diffusion is spherically symmetric. This yields the 
follo wing equation for the energy distrib ution of the cosmic 
rays (JAtovan et al.lll995l : IProfumdlJOOsI ) , 



9/ 
dt 



Di-f) d ^,df ^ 9(P(7)/) , , 



r^ dr dr 



dj 



(1) 



where D{'y) is the diffusion coefficient, P(7) is the energy 
loss rate of the particles, and Q is the source term. Ion- 
ization, bremsstrahlung, inverse Compton, and synchrotron 
radiation draw energy from the cosmic rays as they travel, 
but for cosmic rays younger than 10^ yr the final two pro- 
cesses dominate. In this regime we have 



P(7) = p27^ ~ 5.2 X 10"^S^s"^ 



(2) 



and a simple solution to Eq. \T\ follows for an instantaneous 
burst of cosmic rays at the origin. 



/(^■.*.7) = ^^(1-P2t7r-^ f— )"e-'^/-'>^ 



(3) 



for 7 < 7cut = {P2t)^ (otherwise / vanishes). In this equa- 
tion A''o is the total number of cosmic rays produced with 
a spectrum of dN/d'y oc 7~"; here we take a = 2.2. The 
diffusion length is given approximately by 
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Table 1. Parameters of the diffusion calculations for the two 
scenarios. Dio is the value of the diffusion coefficient at 10 GeV, 
and Di is value at 1 GeV if one extrapolates po wer-law behavio r 
at high energies to 1 GeV (for comparison with |Profumdl2008r . 
Both are in units of cm^/s. The pair-production efficiencies for 
pulsars is one-tenth that of magnetars in both cases. 
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Figure 1. The spectrum of cosmic rays from a single pulsar 
(green dashed curve) , magnetar (red solid curve) or instantaneous 
burst (black dotted curve) at a distance of 200 pc and an age of 
20,000 years under the "LOW" scenario. The total energy release 
is the same for the three models. 
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In contrast with IProfumd (|2008l ). we do not choose the dif- 
fusion coefficient D{'y) to be a strict power la w but assume 
rathe r that it saturates below several GeV (JAtovan et al.l 
Il99d ) 



D(7) = Do (1 + 7/7.) 



(5) 



where we take 7* — (3GeV)/meC^. In addition to giving the 
value of Do, Table[T]presents the values of ^(7) at 10 GeV as 
Dio and the value at 1 GeV if one extrapolates the power- 
law behavior at high e nergy downward. The latter is for 
comparison with IProfumg (|2008i" ). who uses a strict power- 
law relation. 



2.3 Distribution 



We adopt the analysis of iGiU fc Hey J (|2007l ') to determine 
the expected spatial locations of magnetars and pulsars. 
This distribution has two components. One owes to th e 
high mass stars throughout the Galaxy (|Reedll2000l . l2005h . 
This distribution decays exponentially from the center of the 
Galaxy with a scale length of 2.8 kpc and vertically away 
from the galactic disk with a scale height of 45 pc. The 
second component is the overabundance of high mass stars 
known as the Gould belt. Many of the nearby middle-aged 
neutron st ars are probably as sociated with this recent star 
formation (jPopov et al.ll2003l ). Including this population is 
critical because nearby young sources can dominate the ob- 
served spectrum of cosmic rays. To determine the ages of the 



sources, we assume that the birthrate of pulsars and mag- 
netars has been uniform in time throughout the Galaxy, so 
we simply choose 20,000 ages from a uniform distribution 
between zero and 10® years for pulsars and 10^ years for 
magnetars. 

Although we examine sources with ages up to 10^ years, 
the bulk of the emission emerges within the first 10,000 years 
for magnetars or shorter for pulsars. Therefore, even with a 
kick of up to 2,000 km/s, the sources move at most twenty 
parsecs from their birthplaces during their active phases. It 
is reasonable to neglect the motion of the sources through 
the Gala x y dur ing their active period. This contrasts with 
IProfumd ( 20081 ) who used present day positions of pulsars, 
which leads to a more diffuse distribution (e.g. more ex- 
tended perpendicular to the Galactic disk) than when the 
sources were active because of natal kicks given to neutron 
stars. 



3 RESULTS 

There are two important questions to address. The first is 
whether the magnetic fields of the magnetars carry suffi- 
cient energy to account for the observed fiuxes of cosmic 
rays. The second is whether or not magnetars can account 
for this observed spectrum of cosmic rays incident upon the 
Earth, given the assumed spectrum of cosmic ray produc- 
tion. We are interested in both the mean cosmic ray abun- 
dance over the realizations but also the range. Understand- 
ing the range is crucial because one or two nearby young 
sources can dominate the abundances, and we may simply 
be living in a somewhat special place at a somewhat special 
time. It is important to mention that our r esults even for 
pulsars do not agree with those of IProfumd (|2008l ) because 
of the aforementioned differences in modeling. 

Fig. [2]shows the expected abundance of cosmic rays for 
pulsars and magnetars for the two different diffusion scenar- 
ios. Although the total energy reservoir is the same for a 
particular pulsar or magnetar, the magnetars are ten-times 
less numerous, so they are assigned a ten times larger effi- 
ciency so that the total cosmic-ray production is the same 
for the two models. The efficiency is set for the two diffu- 
sion scenarios by insisting that the model predictions agree 
with the Fermi data around 300 GeV (see Tab.[T}. This also 
where the pulsar and magnetar models cross. 

For both diffusion scenarios and both populations the 
resulting cosmic ray abundance comes from nearby old ob- 
jects and more distant young ones. The Galactic distribution 
of the sources goes from approximately three-dimensional to 
disk-like around a distance of 1 kpc. Given the abundance 
of pulsars and magnetars, the youngest source within this 
distance is typically 10* yr or 10^ yr for pulsars or mag- 
netars respectively. This yields cutoff energies of approxi- 
mately 31 TeV and 3.1 TeV, roughly where the flux begins 
to decrease in the "LOW" scenario. Of course, the cutoff 
energy does not depend on the diffusion scenario. In the 
"LOW" scenario the diffusion length is about 1 kpc near 
the cutoff so there are few young sources to contribute to 
the flux at Earth at these energies, and the flux drops and 
recovers at high energies as the diffusion length increases. 

On the other hand, in the "HIGH" scenario the diffusion 
length is a few kpc, so we receive cosmic rays from further in 
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Figure 2. Spectrum of cosmic rays from pulsars (green dashed 
curves) and magnetars (red solid curves) within the Galaxy un- 
der the "LOW" and "HIGH" diffusion scenarios. For each type 
of object, the shaded region gives one-standard deviation vari- 
ation over the realisations. F urthermore, the bold curve gives 
the mean plus the results of iDelahave et al.l 1 120101) excluding 
the contribu tion of pulsars. T he solid triangles give the results 
from A TIC ijChan ^ et al.||2008l). the circles give the results from 
HESS IIAharpnian et al. |l200q) and the squares give the results 
from Fermi JAbdo et al.ll2009f) . For each scenario the pair pro- 
duction efficiency is set so that the total cosmic-ray abundance 
does not exceed that observed by Fermi around 300 GeV. 



the disk from younger sources even around the local cutoff 
— as the older, nearby sources cut out, the younger more 
distant sources seamlessly take over. The transition between 
these two populations leaves an imprint on the expected 
spectrum. The energy of this imprint depends on the details 
of the diffusion process and the age of the sources. As the 
diffusion becomes less rapid, and the objects become older, 
the transition energy decreases. For the magnetar popula- 
tion in the "LOW" scenario, this transitio n occurs tantaliz- 
ingly close to the dip in the HESS data. lAharonian et al.l 
([20091) interpret this dip as a large steepening in the power- 
law around 1 TeV consistent with the magnetar results. The 
pulsars do not follow the dip as sharply but they better fit 
the final HESS data point at about 5 TeV. 

Fig. O shows the expected anti-matter fraction in 
electron-positron cosmic rays for the two different diffu- 
sion scenarios compared and combined with a model that 
excludes a contributi on from either pulsars or magnetars 
(jDelahave et aLluOlOT ). The contribution of pulsars or mag- 
netars increases the expected positron fraction at all ener- 
gies. The effect of magnetars is stronger at lower energies 
than pulsars because the positrons received from magnetars 
are typically older and less energetic than those from pul- 
sars. The variation in the expected positron fraction over 
realisations is larger for the magnetars because the observed 
cosmic rays are produced by fewer dominant sources, so it 
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Figure 3. The observed positron fraction of cosmic rays from 
pulsars (green dashed curves) or magnetars (red solid curves) 
combined with a standard model for cos mic-ray production and 
propagation exclusive of neutron stars from lDelahave et al.l 1120101 ) 
(given by the black curves) under the "LOW" and "HIGH" dif- 
fusion scenarios. For type of object, the shaded region gives 
one-standard deviation variation over t he realisations. The solid 
squares give the results from PAMELA IIAdriani et al.ll2009l) . 



is problematic to access which model does better with a sta- 
tistical tests. At higher energies the pulsars predict a larger 
positron fraction than magnetars, although here the frac- 
tions agree within the expected variation. The disagreement 
with the observations below 10 GeV although statistically 
significant is probably not important because of solar mod- 
ulation of the positron and electron fiux at low energies that 
has not been modelled here. 



4 CONCLUSIONS 

We have derived the expected spectrum and composition of 
cosmic rays from magnetars under the reasonable assump- 
tions that the decay of the external magnetic field powers the 
cosmic-ray production and that the birthrate of magnetars 
is about ten percent of that of pulsars. We assume that pul- 
sar spin-down powers the cosmic-ray production from these 
more weakly magnetized neutron stars. We find that either 
the ensemble of magnetars or that of pulsars can reproduce 
many of the puzzling features of the spectrum of electron- 
positron cosmic rays recently discovered by ATIC and HESS. 
The differences in the appearance and location of these fea- 
tures in the cosmic-ray abundance between the models re- 
sults mainly from the the relative rarity of magnetars. Al- 
though the relative heights of the two bumps in cosmic rays 
detected at Earth do depend on the the assumed power-law 
production spectrum, the location of the transition between 
them does not. On the other hand, the ensemble of pulsars 
cannot easily reproduce the sharp dip in abundance around 
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2 TeV found by HESS - pulsars do reproduce the upturn at 
5 TeV, so a combination of the two ensembles may explain 
both the dip and the upturn. Both pulsars and magnetars 
can account for the positron abundance in cosmic rays with 
energies of about one to one-hundred GeV. With additional 
modelling of the cosmic-ray diffusion, energy losses and pro- 
duction, even more detailed agreement with the observations 
can be achieved. 
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